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Abstract 

We propose a holographic model of baryon interactions based on non- 
supersymmetric D7— D7 flavor branes embedded in the Klebanov-Strassler back- 
ground. The baryons are D3-branes wrapping the 5^ of the conifold with M 
strings connecting the D3 and the flavor branes. Depending on the location of 
the latter there are two possibilities: the D3 either remains separate from the 
flavor branes or dissolves in them and becomes a flavor instanton. The leading 
order interaction between the baryons is a competition between the attraction 
and the repulsion due to the a and uj mesons. The lightest O^"*" particle o" is a 
pseudo-Goldstone boson associated with the spontaneous breaking of scale in- 
variance. In a certain range of parameters it is parametrically lighter than any 
other massive state. As a result at large distances baryons attract each other. 
At short distances the potential admits a repulsive core due to an exchange of 
the Lo vector meson. We discuss baryon coupling to glueballs, massive mesons 
and pions and point out the condition for the model to have a small binding 
energy. 



1 Introduction and overview 



Understanding nuclear forces is an important problem with many interesting questions 
still not answered. One of them is the smallness of the nuclear binding energy. Despite 
the strength of the strong interactions the force between two nucleons is rather weak: 
the binding energy is only a few percent of the typical baryon mass. At long distance 
the force is mediated by exchanges of the lightest isoscalar particles:^ a (0"''"'") and 
Lo (1 ). The masses of a and uo are close, hence they effectively cancel each other's 
contribution to the binding energy. The origin of such a coincidence is not easy to 
establish. Moreover, the identification of the scalar meson and its relation to a pion- 
pion resonance is a long-standing puzzle. Theoretically one may also worry about 
the glueballs though they are heavier than a and in many cases can be neglected. 
Accounting for these phenomena using first principles of QCD has been found to be 
a non-tractable task. 

Provided the approximate equality of a and uj masses is not totally coincidental, 
to find a reason it would help to study a toy model, in which the contributing factors 
can be considered separately. A natural idea would be to look for a holographic model 
of interacting baryons, in particular because in the large N^. limit, required by the 
holographic approach, the exchange by many particles as well as by glueballs will 
be suppressed. This would leave to understand whether the contribution from the 
single scalar and vector meson exchanges approximately cancel each other. Attempts 
have been made recently to find such a holographic model. In [1, 2, 3] the baryons 
were considered in the context of the Sakai-Sugimoto model [4]. It turns out that 
in this model the attractive force between the baryons is always weaker than the 
repulsive one [3]. The light scalar meson happens to be heavier than the lightest 
vector meson resulting in a net repulsive force. Clearly this is very different from 
hadronic physics. In order to shed some light on the nucleon binding energy and 
in particular its smallness, we should first try to find a model featuring a binding 
potential and then try to examine whether it is naturally small. 

Motivated by the above considerations, in this paper we study a holographic 
model based on D7— D7 branes of [5] embedded in the Klebanov-Strassler (KS) back- 
ground [6]. This system is scale invariant in the UV. Adding the D7— D7 branes 
spontaneously break conformal symmetry by a vev of a marginal operator at some 

^We implicitly assume the case of baryonic matter when the contribution due to one-pion exchange 
is believed to average to zero. 
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scale tq. When this scale is larger than the internal scale of the gauge theory = e^^^, 
the lightest scalar meson is parametrically light as a pseudo-Goldstone boson of the 
conformal symmetry. This meson gives the leading contribution to the attractive 
force and we will retain the notation a for it in the remainder of the paper. Similarly, 
we will call u the lightest vector meson, which corresponds to the fluctuations of the 
f/(l) gauge field on the world-volume of the D7— D7 branes. 

The model in question has the following hierarchy of light particles. We are 
working in the quenched approximation, that is the number of D7— D7 branes Nf is 
much smaller than the number of colors M. Hence we ignore the backreaction of the 
open string sector. In particular, the mass of glueballs remains the same as in the KS 
theory and therefore is ro-independent. The typical scale of the glueball mass is 

mgb~^, X = gsM. (1.1) 

In the regime Tq ^ the theory is (almost) conformal and therefore the mass of 
mesons can depend only on the scale of symmetry breaking tq 

a'X 

The pseudo-Goldstone boson cr is parametrically lighter 

m<^~-|mgb. (1.3) 

^0 

We will explain this result later in the text. 

As ro approaches the mesons become lighter, while the pseudo-Goldstone grows 
heavier. Around the minimal value ro = r^ all mesons have approximately the same 
mass of order rrigb. This is the most interesting regime of parameters because for 
ro ~ r^ the approximate cancelation of the attractive and the repulsive force can 
occur naturally. Hence it is crucial to know the precise hierarchy of the lightest 
particles in this limit. Such a detailed computation is beyond the scope of the present 
work. A recent publication [7] based on the approximate analysis ignoring an intricate 
dynamics of the gauge fields on D7— D7 branes claims the following spectrum in the 
limit ro = r^ 

< mQ++ < < mi++ . (1.4) 

Here 0"'""'" and l"*""*" denote the lightest glueballs with the corresponding quantum 
numbers and we have used the results on the KS glueballs from [8, 9]. This result is 
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encouraging and we will investigate the phenomenology of the model assuming that 
< for tq ~ r^. 

A baryon in our setup is represented by a D3-brane wrapping the of the conifold 
and a set of M strings connecting it to the D7— D7 branes. For large values of the 
parameter tq ^ the string tension is smaller than the force exerted on D3 due to 
curved geometry. To minimize the energy D3-brane will settle near the tip of the 
conifold at r ~ with the D3— D7 strings stretched all the way between r and tq. 
When ro is significantly close to the geometry can be effectively approximated by 
the flat one and creates only a mild force. The string tension wins, and the D3-brane 
is pulled towards the D7— D7 branes and dissolves there becoming an instanton. 

The description of the baryons in the two regimes is essentially different. Never- 
theless, in both regimes the D3-brane is close to the tip of the conifold and the form 
of the glueball induced potential is approximately the same 

-mgb|a:| 

\/,b~±— ^. (1.5) 

Here x is the 3-dimensional distance between the two baryons which is assumed to 
be much larger than the size of the baryons p. In the above calculation we estimated 
the force at large distances, where it is dominated by the lightest particles. In this 
approximation baryon is treated as a localized source. The sign in (1.5) depends 
on the spin of the exchanged particle: the scalar induces attraction while the vector 
yields repulsion. In the KS theory the lightest O"*"^ glueball is lighter than the vector 
l^"*" and the coupling constant can in principle be read off the explicit form of the 
wavefunction [8, 9]. Although the glueballs will generate some potential, this fact will 
not be important for us, since the glueball potential is suppressed by Qs as compared 
to the potential generated by mesons. 

Let us first focus on the regime ^ r^, in which one adopts the string picture 
of the baryon. The main interaction at large distances is due to the lightest mode a. 
The effective potential of the scalar meson exchange reads 

^scalar ~ \ 7 " (1-6) 

5fslog(ro/r,) \x\ 

It is indeed parametrically larger than (1.5) in the holographic regime Qs — )■ 0. More 
precisely, to ignore glueballs we need to satisfy gs\og{rQ/r^) <^ 1. 

Since the cu-meson is much heavier than a in the r ^ regime, the repulsive 
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interaction is subdominant at large distances, although the potential is similar to (1.6) 

I Q-m^\x\ 

5(slog(ro/rJ |x| 

Theoretically it is possible that in this regime the potential will stay attractive 
even at distances |x| <C i^Z^- Phenomenologically this is not an appealing scenario. 
Assuming that the repulsion will start winning at some |a;| ~ m^^ this will imply the 
binding energy of the order 

-E'binding ~ — f^vcctoT ~ r, " , • (1-8) 

9 s g^Ma' 
This is parametrically smaller than the baryon mass, 

but the relative smallness is due to the artifact of the holographic approach - the 
large 't Hooft constant A = gsM ^ 1, not a cancelation between the attraction and 
repulsion. 

Now let us consider the opposite regime tq — )■ r^. For tq small enough the wrapped 
D3-brane will dissolve in the D7— D7 and will be described by an instanton. We will 
assume that the scalar a is still lighter than the vector to in this regime. When 
ro reaches its minimal value r^, the resulting configuration of the D7— D7 branes is 
invariant under an emergent symmetry. The wavefunction of a is odd under 

1j2 G f/(l) and therefore the leading coupling of a to baryons vanishes. The same 
phenomena also occurs in the Sakai-Sugimoto model [3] . By varying ro near the point 
ro = r^ one can tune the coupling of a to be small. 

The net potential in this case can be written in the form (cf. [10]) 

'^ = ^U-T^-sl-TT-]- (1-10) 

It is valid only for |x| large enough. Here ga and g^^ encode the coupling of the mesons 
to the baryons and we have kept the factor of ga explicit. If m„ < m^, the potential 
is attractive at large distances no matter what the couplings are. On the other hand 
if g^j is small enough, at distances shorter than the vector interaction "wins" and 
the potential becomes repulsive. The binding energy 

r 

^binding ~ «:^PJ^ , (1.11) 
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is suppressed by a small dimensionless number k, which is related to the smallness of 
the coupling g„ and the fact that mo- and are of the same order. This suppression 
is in addition to the relative factor (compared with the baryon mass). The extra 
factor K is phenomenologically promising as it represents the near-cancelation of the 
attractive and repulsive forces responsible for the small binding energy in hadron 
physics. 

This paper is organized as follows. In section 2 we describe the setup and de- 
termine the mass scales of the mesons and glueballs. In section 3 we introduce a 
baryon into the D3— D7 geometry and discuss its equilibrium positions. There are 
two limiting cases of interest: for large tq the baryon sits at the tip of the conifold, 
while for r ~ e it dissolves into D7— D7 and becomes an instanton. In section 4 we 
investigate the relevance of the glueballs to the interaction of the two baryons. We 
find that it is suppressed by Qs compared with the interaction carried by the mesons, 
which we study in section 5. We consider two regimes, corresponding to a stringy 
and an instanton pictures of baryons. Our analysis is summarized in section 6, where 
we also discuss the implications of the model and open questions. In the appendices 
we review the question of the stability of the D7— D7 embedding (appendix A) and 
obtain estimates for the mass and the size of the baryon (appendix B). 

2 Mesons 
2.1 The model 

Let us start with the description of the holographic setup. To simplify the presen- 
tation we will frequently compare our model with the more familiar one of Sakai 
and Sugimoto [4]. The KS geometry results from D3-branes and M fractional 
D3-branes placed on the conical singularity of a conifold, a cone over the T^'^ space. ^ 
The dual gauge theory is a A/" = 1 SYM with the gauge group SU{N + M) x SU (N) 
and bifundamental matter in a left and right doublets of the global SU{2) x SU{2) 
symmetry. The theory flows to the IR by undergoing a chain of Seiberg dualities and 
ends as a confining SU{M) SYM. The light states in the low-energy theory are super- 
glueballs and hybrids containing bifundamentals. They were extensively studied in 
the SU{2) X SU{2) singlet sector [8, 9, 11, 12] and beyond [13, 14]. 

The Sakai-Sugimoto model is based on the A^^ D8 and D8 flavor branes immersed 

^T^'-^ ~ 5^ X 5^ and fractional D3-branes are the D5-branes wrapped on S'^ of T^'^. 
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in the geometry of D4-branes compactified on a circle with anti-periodic boundary 
conditions [15]. The D8 and D8-branes smoothly connect to form U-shape configu- 
rations thus realizing the dynamical chiral symmetry breaking U{Nf) x U{Nf) — )■ 
U{Nf). Similar U-shape configurations of D7— D7 branes in the context of the coni- 
fold were first proposed in [5] for the Klebanov-Witten (KW) geometry [16], a con- 
formal M — limit of the KS model. They were extended to the KS case in [17]. 
Now we will review this construction in more detail. 

The original and fractional D3-branes are placed at the tip of the conifold spanning 
the Minkowski directions. By analogy with the Sakai-Sugimoto model four of the eight 
directions of the D7 and D7 are in Minkowski, three are wrapping the of the T^'^ 
and the remaining one is a line in x S*^. The geometry of the profile is controlled 
by the effective DBI action 



where is the induced metric and J-" is the U{Nf) gauge field on the world-volume of 
the D7 and D7. The dilaton is constant e*^ = gg, unlike in the Sakai-Sugimoto model. 

In both the KW case [5] and the KS case [17] the position of the D7 and D7 on the 

X S"^ can be restricted to the equator of S"^ by an appropriate SU (2) symmetry 
transformation. Effectively one obtains a U-shaped profile of the connected D7 and 
D7 pairs on a cigar-like manifold R'^ x S^, see figure 1. 

To make a better connection with the Sakai-Sugimoto construction we map the 
parameter of the two models as follows. The M^k mass scale of the Sakai-Sugimoto 
model, which is related to the coordinate distance Ukk to the cigar from the origin, is 
mapped to the KS scale = e^^^. The latter is, roughly speaking, a similar distance 
measuring the deviation of the deformed conifold from the singular one. The minimal 
distance Uq between the D8— D8 branes and the and the origin U = is mapped to 
the minimal distance from the origin to the D7— D7 branes vq. The case tq = is 
the antipodal configuration, for which the position on the of the cigar (equator of 
the 5*^) is the same for all r. 

In the KW case the regularized energy of the D7— D7 configuration does not 
depend on tq. Therefore there is a massless particle corresponding to variations of r^. 
It is a Goldstone boson of the spontaneously broken scaling symmetry. The symmetry 
is broken dynamically in the IR, rather than by the boundary conditions in the UV. 
Indeed the asymptotic angular separation along the equator A0(roo) = -\/67r/4 is the 
same for all tq. The angle is dimensionless and the boundary conditions at infinity 




(2.1) 
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(a) (b) 

Figure 1: Embedding of tiie D7-branes. The branes wrap and span a line in R'^ x S^. 
The lowest point of the D7— D7 profile is at tq. In the antipodal case (a) the configuration 
descends all the way down to the tip of the conifold ro = r^. It is smooth at the tip. In the 
non-antipodal case (b) the lowest point of D7— D7 is at some distance from the tip, vq > r^- 



(j) — const do not break conformal symmetry. This is in contrast with the Kuperstein 
embedding [18] (the KS analog of the Karch-Katz embedding [19]) specified by the 
dimensionful parameter 2:4 = /i. 

In the KS case the conformal symmetry is also broken explicitly by the geometry. 
Therefore the variation of tq is no longer massless and becomes a pseudo-Goldstone 
boson. As we demonstrate below, the pseudo-Goldstone boson can be made signifi- 
cantly lighter than other mesons and glueballs. 

Before proceeding to the discussion on the particle spectra in the next section we 
will briefly discuss the unbroken symmetries of the model. The original KS geometry 
has an SU{2) x SU{2) symmetry that rotates T^'^. The left SU{2) is preserved by 
D7— D7 as it rotates the 5*^ of the conifold, which the D7— D7 branes wrap. The right 
SU{2) rotates the S"^ and is therefore broken. In the antipodal case the U{1) of the 
second SU{2) survives with important consequences for the phenomenology of the 
model. 

Another symmetry of the KS solution is the X-symmetry: a Z2 which interchanges 
the two S"^ of the conifold with simultaneous flipping of the F3 and signs [16]. In 
the dual field theory this symmetry acts as charge conjugation combined with the 



interchanging of the bifundamentals. The D7— D7 branes will break this symmetry, 
because they wrap a particular S^. For a generic r, the is not symmetric with 
respect to the swap of the two-spheres, except for the point at the tip of the conifold. 

The D7— D7 branes break supersymmetry and may significantly change the KS 
model through a backreaction. In this paper we assume the quenched approxima- 
tion Nf <^ M when the number of flavor branes is small. Thus we can ignore the 
backreaction of the D7— D7, in which case the glueball sector remains approximately 
supersymmetric, while the flavor sector has broken supersymmetry. To simplify the 
presentation we sometimes will omit the Nf factor in what follows, since the depen- 
dence on Nf can be easily restored. 

2.2 Masses of vector and scalar mesons 

In this section we will estimate the masses of lightest vector and scalar mesons in 
the models of [5, 17] and show that the pseudo-Goldstone boson a can be made 
parametrically lighter than the vector u. This is in a stark contrast with the Sakai- 
Sugimoto model, in which a is always heavier than u [20]. To see how this works 
we start with the conformal background of the KW solution and embed the D7— D7 
branes as described in the previous section. The embedding can be parameterized by 
4>{r), namely, the functional dependence of the angle along the equator of S"^ on the 
radius. The effective action for 0(r) has the form [5] 

S^-L^jr^f^^^. (2.2) 

The explicit general solution can be readily obtained. It is U-shaped and parame- 
terized by the minimal radius tq of the profile. For any tq the asymptotic angular 
separation between the branches is the same A0 = ^vr. To understand why this 
is the case, we need to remember that the KW theory is conformal. The D7— D7 
configuration breaks the conformal symmetry spontaneously by introducing the only 
dimensionful parameter tq, while the boundary conditions at infinity (p ~^ constant 
do not break the symmetry. As a result the mode d(j)/dro is a massless Goldstone 
boson that does not change the boundary conditions. Besides the Goldstone boson of 
the conformal symmetry (and the pions - the Goldstone bosons of the spontaneously 
broken chiral symmetry) all other meson modes, including the vector u, acquire a 
mass of order 
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where is the glueball mass scale 

g2/3 ^ 

m,, = — = (2.4) 

The combination e^^^ /a' is identified with the strong couphng scale A. In the regime 
tq ^ the mesons are heavier than the glueballs. 

In the KS case the conformal symmetry is broken explicitly by the deformation 
parameter e. In this case rg is still a free parameter tq > r^, but the angular separation 
depends on tq. It is maximal A0 = tt for tq = and approaches the KW value 
A(/) = for To ^ [17]. In the regime tq 3> we can treat the effect of as 
a small correction. Hence the masses of regular mesons remain essentially the same, 
specified by the scale tq, but the massless scalar mode a will acquire a small mass of 
order 

^2/3 ^4/3 ^2 

m„ ~ TT = -t: m„b • (2.5) 

We will explain this result momentarily. A priori it is not obvious why the mass of a 
should be positive, i.e. why the resulting configuration should be stable. To under- 
stand why this is indeed the case we would like to further investigate the connection 
between the asymptotic angle A0 and tq. Let us for a moment return back to the 
conformal case = and consider a probe configuration of the following form: a 
straight line along the radial direction r from infinity until tq at some fixed (p = (pi, 
an arc at fixed r = Tq between 0i and 02 = 0i + A0 and a straight line at fixed 
= 02 from To to infinity. The action of such a configuration as given by (2.2) reads 

^N = ^ + aro^(A0-A0o). (2.6) 

Here a = 1/ {G^^'^gsCt'^) is some positive constant and A0o = "\/6/2. Let us analyze this 
result. The effective action (2.6) is not bounded from below if A0 = 02 — 0i < A0o. 
In this case the brane configuration is unstable and the U-shaped stack is pulled "up" 
towards large r and disappears. The D7 and D7-branes annihilate. When A0 > A0o 
the potential is bounded from below and the stack stretches down to minimize tq. The 
resulting configuration is stable. Finally, if A0 = A0o the parameter tq is marginal 
and can acquire any value. In our analysis we ignored the fact that (2.6) is based on 
some probe configuration, which does not minimize the action even locally. Hence 
what we perceived as a stable configuration, minimizing tq, may not in fact be such. 
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From an explicit calculation we infer that the correct angular separation in the KW 
case is A0o = "\/67r/4. This value is larger than the approximate one (-\/6/2), which 
is consistent with the logic that (2.6) accurately predicts the fate of the unstable 
configurations but may not be accurate for the stable ones. 

The effective action (2.6) has a simple field theory interpretation. The dimension- 
less parameter A0 — A0o is a coupling constant in front of a marginal operator 
The latter acquires a vev ($^) ~ r^/a"^. Indeed the non-zero value of tq is attributed 
on the field theory side to a vev of some chiral symmetry breaking operator [5]. 

What happens with the effective action when we turn on e? Then we have dif- 
ferent To for different asymptotic angular separations A0. Qualitatively this can be 
described by the following effective action 

5eff [$, A0] = a($2 - $g(A</.))^(A0 - A0o) , A0o = ^vr , (2.7) 

and the vacuum value $o(A0) is related to ro(A0). It is clear that for A0 > A0o the 
configuration of the D7— D7 branes extends until $ = ro(A0)/a', which is a stable 
position with the positive mass of small fluctuation. 

One should not worry about stability of the D7— D7 profile with respect to the 
SU{2) charged fluctuations, i.e. the Kaluza-Klein modes on the S^. This is because 
all such modes are massive in the KW case with the mass specified by rg as in (2.3). 
(See the details in appendix A.) In the tq ^ e^/^ regime these modes will of course 
remain massive. 

The argument above could be made more rigorous if we reformulate the question 
of the sign of as a question of stability of the classical configurations found in [17]. 
It was explained there that the D7 action is minimized when the world- volume gauge 
field vanishes. As a result we can discuss stability of the U-shaped D7— D7 configu- 
ration in the geometry of the deformed conifold with no flux present. As was shown 
there the classical solution with given boundary conditions (asymptotic angular sep- 
aration A0) is unique and hence must be stable. Indeed there is no singularity at 
small r for the D7-branes to end there. Also, the analysis in the KW case shows 
that for A(j) > A(f)o the D7 and D7-branes do not annihilate. Hence for the given 
angular separation the unique solution found in [17] must be the minimum of the 
action, not just an extremum. Let us notice here that a similar argument does not 
work for the U-shaped configuration of the Kuperstein D7— D7 brane embedding [18] 
(or Karch-Katz D7— D7 configuration in the AdS^ x space [19]). Although in that 
case there is a unique classical solution for the separation at infinity, the asymptotic 
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angular separation is zero (the branes are parallel to each other at large r) and the 
configuration is pulled to infinite tq. 

To calculate the mass of a in terms of tq we need to know A0 as a function of 
ro and accompany the potential (2.6) by a kinetic term. To estimate A0 we expand 
the corresponding expression (it can be found in [17]) at large radius. The leading 
deviation from the KW geometry is given by the deviation of 

K{t) sinh r = (sinh r cosh r — r)^^^ 

from r^/e^/^. Expanding we get 

i^^(r)sinhr^-^(^l-— log^j. (2.8) 

As a result the leading deviation of A0 from A0o is 

A0 = A0o-&4log^o, (2.9) 

where h is some unimportant numerical coefficient. 

Now we need to calculate the kinetic term. We substitute Sr^ for the field $, 
since the geometrical meaning of the latter is the variation of ro by (5ro. To evaluate 
the kinetic term we simply plug the time dependent wavefunction </)(r, t) into the 
effective action and expand to quadratic order around the classical solution for 0. In 
the leading order one can restrict to the zero mode solution of the conformal case, 
i.e. the variation d/dvQ of the classical solution 

S<l>{r) = -5ro{xo)^\L^,. (2.10) 
ro y r^ — Tq 

Hence the full action, up to unspecified numerical coefficients A and B, reads 

'-^I'"' (^toMa')^log.„ (^)' - B J los^o (fro)'^) , (2.11) 

and the mass of pseudo-Goldstone meson is given by (2.5). 

This result has a clear interpretation in terms of an effective field theory. In the 
regime when ro ^ r^ we have a separation of scales. At energies much smaller than 
Muv ~ $ = ro/a' all massive degrees of freedom associated with the D7-brane can be 
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integrated out, except for the light pseudo-Goldstone boson. The low-energy effective 
field theory is the original KS theory perturbed by some operators 



There is also a kinetic term for the pseudo-Goldstone boson (9$)^. Comparing (2.12) 
with (2.11), we conclude that the expansion starts with the operator of dimension 
6, i.e. it is irrelevant. The absence of relevant perturbations can be attributed to the 
quenched approximation, which assumes that the field theory returns to the original 
one once the D7-brane is completely removed {tq — > oo). This is in contrast with the 
pseudo-Goldstone boson of chiral symmetry in QCD, where the chiral symmetry is 
broken explicitly by a small mass nig ^ Muv = ^qcd 



Here the operator \E'\& of dimension 3 is relevant and therefore the effect of (2.13) on 
the low-energy physics increases with Aqcd — ^ oo. 

3 Baryons 

Following the general prescription of Witten [21] we introduce a baryonic vertex (b.v.) 
as a D3-brane wrapping the cycle of T^'^. Because of the flux of the 3-form 



D3-brane carries M units of charge under the U{1) field on the brane. This implies 
that there must be M strings ending on the wrapped D3 brane. The other ends of 
the strings are attached to the flavor D7— D7 branes, representing M quarks. In this 
way the baryon gets flavored. The original setup [21] was in AdS^ x 5*^ which is 
conformal non-confining background and the strings ended on the boundary. This 
is a holographic representation of a baryon made of infinitely massive quarks. A 
similar discussion for confining background was presented in [22, 23]. Dynamical 
baryons, as opposed to infinitely massive ones, were introduced by ending the strings 
on flavor branes rather than on the boundary. This was done in the context of the 
Sakai-Sugimoto model in [24, 25]. 

The strings pull the b.v. towards the D7— D7 branes. This is opposed by the 
geometry which, through the DBl action, attracts the b.v. to the tip of the conifold. 




(2.12) 



6C = mg^^ . 



(2.13) 




(3.1) 
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In the case of the Sakai Sugimoto model the strings always win [25]. As we are going 
to demonstrate, this is not the case for the KS background. 

The energy of the b.v. at a point r = e^^'^ cosh^^^ r is given by (here we use the 
notations of [6, 26]) 

[ e-*v/=^~^/ii/2(^)coshr. (3.2) 
a' Js-i QsOi' 

Here /i(r) is the warp-factor of the KS geometry. For large r the energy scales as 

M r 

^b.v. ~ — rlog-. 3.3) 

a' 

The energy of the M strings is given by 

1 ^ M n e^/^ 

a' 

which for large r behaves as 



i?stnngs = M-/ v/^^dr~-/ — -^dr, (3.4) 



M 

^strings 7 (^0 - ^) • (3.5) 

a' 

The position of the D3 along the radial direction will be determined by minimizing 
the total energy 

M 

^tot = — (^rlogr + 5|ro-r|)), (3.6) 
a' 

where A and B are some numerical coefficients. Hence the b.v. will sit at r ~ 
riR. Keeping in mind that there are M strings which can affect each other through 
backreaction, it is hard to find the actual values of A and B. Nevertheless, it is clear 
that for large tq the D3 will prefer to settle in the IR and the baryon will correspond 
to M strings stretching between the D3 and D7— D7 (figure 2). 

In such a case the properties of the baryon are determined by a composite system 
of the b.v., the M strings, their endpoints on the b.v., and their fiavored endpoints 
on the D7— D7 branes. The endpoints on the b.v. are fermionic thus assuring the 
antisymmetric product of M fundamentals [21]. On the other hand the fiavored 
endpoints have to be bosons so that one can recast using the Hartree approximation 
the well known picture of large M baryons [27] . In this picture the size of the baryon 
scales with and the mass as M. 
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Figure 2: (a) Baryon as a D3-brane wrapped on the S*^ of the conifold. It sits between the 
tip of the conifold (r^) and the lowest point of the D7— D7 configuration (ro). M strings, 
stretched between the D3 and the D7, pull the D3 up, but the geometry force is stronger 
and keeps the D3 close to the tip of the conifold. (b) If the D7— D7 configuration descends 
close enough to the tip, the string tension wins over the gravity and the D3 dissolves in the 
D7 becoming a fiavor instanton. 



For small tq the space is locally flat, the energy of the wrapped D3 will be ap- 
proximately independent of the location r and the strings can win. Expanding the 
above integrals around r = 0, one finds 

i^b.v-^ 62/^/^/^(0) + 0(r^)), (3.7) 
a' 

where the warp factor was traded for a dimensionless function h oc {gs^Mc^f' e'^l"^ I{t) ^ 
and 

i^strings ~ ^ e^'^ (|ro - r| + 0[tT)) . (3.8) 

In the region of sufficiently small ro, in which one can approximate the energy of D3 
by a constant, D3 will be pulled towards the D7— D7 and dissolve in them, i.e. it will 
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be equivalent to an instanton configuration of tlie gauge fields on the D7— D7 world- 
volume [28]. In the latter case the instanton coupling comes from the Chern-Simons 
action of the D7— D7 branes: 

where T = dA + ^ A ^ and 

u^iA) = Tr (^-A AJ^^ + ^A^ AJ'+^A'^y (3.10) 
Expanding the U{Nf) field into the SU{Nf) and the U{1) parts 

one obtains from the first term in (3.10) 

Scs-M j d^x dr io Tr F^^ . (3.12) 

Fjnn is the field strength associated with the SU{Nf) gauge fields and the indices 
■m,n = 1,2,3,4 run through the space coordinates and the radial direction. For the 
instanton configuration of the non-abelian field A the result is the instanton number 
n times M units of charge under the U{1) gauge field A. Since the field Aq is dual 
to the baryon number chemical potential in the dual field theory, n is the baryon 
number. 

The effective 5-dimensional action for the gauge fields on the D7— D7 branes con- 
sists of the interaction part (3.9) and a kinetic term, obtained from the expansion of 
the DBI action (2.1): 

5dbi = -^^ / d'^d^M^)Tr[M^)-F^l + 2J-y . (3.13) 
Here a new radial coordinate z is introduced, such that 

dsl = h-'/^dsl, + h'/' {dz^ + r^{z)dnl) (3.14) 

and 

, , , cosh t(z) 

k(z) = 3.15 
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Integration over z in (3.13) is implied over both branches of the D7— D7 branes. The 
solutions to the equations of motion can be found as an expansion in A~^, starting 
from flat-space BPS instantons as a leading order approximation. 

In this approximation the mass of the instanton is given by the expression (3.7): 

mtaryon ~ ^ e'/' = MA . (3.16) 

Higher order corrections come from the non-zero curvature at the position of the 
baryon and the Chern-Simons interaction. A more detailed analysis of corrections is 
done in appendix B. It is very similar to the analysis in the Sakai-Sugimoto setup [29]. 
Here we notice that also in the KS case the size of the baryon remains small, 

g2/3 

P~^, (3.17) 
which makes higher string corrections important. 



4 Glueball induced potential 

Now we turn to the interactions between baryons. First we will estimate the glueball 
contribution to the baryon-baryon interaction potential. To calculate the potential we 
would need to know how glueballs couple to baryons. Notice first that the scalar and 
vector glueballs couple to the baryon in different ways. The former couples through 
the DBI term, while the latter through the Chern-Simons term. Let us start with 
the scalar. We will consider a mode which warps the six-dimensional geometry of the 
conifold 

ds2 = /i"i/2d4i + e^/ii/M52. (4.1) 

Here X is a small perturbation and we are going to find the effective action describing 
its dynamics. In the KS background this mode apparently couples to other glueballs. 
We will neglect this interaction as it will only change the numerical coefficient, but 
not the parameter dependence of the interaction. We start with the action in the 
bulk expanded to the second order in X: 

^buik « ^ / ^^''xbg^^'bRAB + ... (4.2) 

We focus on the singlets with respect to the S\J (2) x SU (2) symmetry as the charged 
glueballs tend to be more massive. Now we use the relation 

bR,, = \ji,,hd''X, (4.3) 
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where d"^ is the Minkowski Laplacian. After integration by parts and evaluating the 
integral along the radial r and angular directions we arrive at the effective action 

^buik[X] ^ (^^) / {{dXf - m\X^^ . (4.4) 

Here we separated the Minkowski and r dependence of the wave function X(x^, r) = 
X{x^^t\}{T\ The four- dimensional mass m| has a typical scale of the glueball mass 
mgb. Its numerical values were computed in [8]. It was also shown there that there 
are modes with '?/^(r) not vanishing at the tip '?/'(0) ~ 1. 

The interaction with the baryon is even easier to find. As we discussed earlier 
the baryon can either be described by a D3-brane wrapping the S'^ cycle or by an 
instanton inside the D7— D7 flavor brane. In any regime it sits close to the tip of the 
conifold r ~ and its mass is given (at the leading order) by Me^^^/a'. The warp 
mode X will change the geometric volume of 5*^ and hence the mass of the baryon, 
leading to the interaction term 

S^^\X\ = / dx°X^(0). (4.5) 

Taking '0(0) ~ 1 we find the glueball induced potential for two baryons separated by 
the distance |x| to be 

-mgb|x| 

V,, ~ . (4.6) 

The derivation of the repulsive vector meson potential proceeds along the same 
lines. Vector glueballs even under charge conjugation were studied in [9]. Contrary 
to the scalar glueballs, vector equations can be diagonalized and one can enjoy using 
the explicit form of the wavefunction. Of particular interest is the lightest vector dual 
to the f/(l)_R current, which appears from the fluctuations of the metric and 5-form 
field. This glueball interacts with the baryon through the coupling of the wrapped 
D3 brane to the four form C4: 



^D3~ — / SC, (4.7) 
a Jd3 

where SC4 can be read off from [9]: 

6C^ ~ QsM^a'^ V Ag,Adg, + ..., g, = + . . . (4.8) 
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Here V = ip{T)V^{x'^) dx'^ denotes the bulk glueball field and the conventional para- 
metrization of T^'^ is used [6]. The fluctuation also affects the metric 

Plugging this into equation (4.2) we arrive at the following expression for the effective 
action for 

^buik[^] ^ g!M' J d'x ((a^K)2 - mlV^^ , (4.10) 

which gives the same interaction potential (4.6) modulo a numerical prefactor. 

As we will infer from the next section, the interaction carried by mesons is en- 
hanced by as compared with the one by glueballs. Therefore we will neglect 
glueballs in the later discussion. Notice however that smallness of Qs is dictated by 
the holographic approach. 

So far we have ignored the massless glueballs of the KS theory [30, 31], which could 
dominate the potential at large distances. There are two exactly massless particles: 
the Goldstone boson of the spontaneously broken U{1)b baryon symmetry and its 
CPT-partner. These glueballs are odd with respect to the X-symmetry, which we 
have discussed above. Since the geometry of the deformed conifold is invariant under 
this symmetry the cycle of minimal volume located at the tip of the conifold is also 
invariant. As a result the baryon wrapping this cycle does not couple to the X-odd 
massless glueballs at the lowest order. The X-symmetry is broken by the presence of 
D7— D7 branes. The D3— D7 strings would pull the D3-brane in a X non-invariant 
location. This in turn would create a coupling between the baryon and massless 
glueballs. The wavefunction of these glueballs behaves as near the tip [31] and 
provided that the D3-brane is sufficiently close to the tip the coupling is not large. 
Although these glueballs will be responsible for the dominant force at large distances 
their contribution to the binding energy will still be Qg suppressed compared with a 
and w mesons. This argument allows us to neglect the effect of these particles making 
our model more similar to hadron physics. 
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5 Meson induced potential 



5.1 Regime of large tq ^ 

In this regime the baryon is the b.v. that sits close to the tip of the conifold and is 
connected with M strings to the D7— D7 branes. The pseudo-Goldstone boson a is 
the hghtest particle except for the massless pions (and rj') and the pair of massless 
glueballs. This fact guarantees the attraction of two baryons at large distances. To 
estimate the binding energy we need to combine the action (2.11) for a with the term 
describing its interaction with the baryon through the D3— D7 strings. 

When ro is large the pseudo-Goldstone boson is simply generated by the dilatation 
d/dro. Hence to find the interaction term we need to estimate how the energy of the 
baryon (and the attached strings) will change when tq is changed. Clearly the main 
contribution comes from the energy of the D3— D7 strings, which depends on tq (3.5). 
Hence the interaction term is 

M 



a' 



dx^Sro. (5.1) 



Together with (2.11) this gives the attractive potential 

I Q-m„\x\ 

^scalar ~ ] fH • (^-2) 

fif^logro |a;| 

As was mentioned before the potential energy is I/qs times enhanced as compared 
with that of the glueballs. 

The vector meson couples to the strings as follows. The endpoints of the strings 
act as sources for the gauge fields on the D7— D7 world-volume. For a static baryon 
we obtain in the point-source approximation 

5'strmg[i] = M j d^x J% ~ M y dx% . (5.3) 

Alternatively this result can be obtained from the boundary term in the action for 
the strings: 



M J A. (5.4) 
The kinetic term for Aq = do{x'^)(f{r) comes from the DBI action on the D7— D7: 
^DBi[a] ~ 9sM^ J d^x {{ddoy - mldl) . (5.5) 
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This gives the same expression (5.2) for the potential modulo the sign and a numerical 
prefactor. 

A more detailed analysis is required in order to understand, whether the repulsion 
will overcome the attraction at short distances. Such a calculation however is beyond 
the scope of the present paper. A more interesting prospect opens in the regime of 
small To, which we discuss in the next section. 



5.2 Regime of small tq ~ 

In this regime the D3-brane dissolves inside the D7— D7 to form an instanton. To 
calculate the meson induced potential at large distances, we need to calculate the 
interaction between the mesons and the instanton solution on the D7-branes. We start 
with the pseudo-Goldstone boson, which in this case has a mass of order mo- ~ "Tigb. 
Its wavefunction can still be parameterized by the deviation of the lowest point of the 
profile 6ro. But in this case there is only one scale tq ~ and it is actually easier to 
estimate the result than in the ^ regime. The wavefunction can be expanded 
into a product of the Minkowski wave function o"(x^) and the internal wavefunction 
ip{4>), which depends on the angular variable along the 5*^. (We choose (p instead 
of the radius because it is a good coordinate near the tip = 0.) The internal 
wavefunction is of the order of ip{0) ~ and a simple estimate gives, cf. (2.11), 

SM ~ ^JfM^ I [^daf - mla^) . (5.6) 

To calculate the interaction term we need to expand the DBI action keeping the 
dependence on the non-abelian gauge fields, 

^na - ^ y d^xd0 rVr^ + 6(rO^Tr [hPf^ + ^^^^,) , (5.7) 

and then expand up to the lowest order in 6r. In the above expression r' = dr/d0. 
Taking into account that the non-abelian fields present a point-like instanton at the 
origin a straightforward calculation gives 

5na ^ / dxV. (5.8) 

Therefore the resulting potential is 

g2 Q-m„\x\ 

Scalar — i i • (5-9) 
A-KQs \X\ 



21 



Here we have introduced a dimensionless coupling constant gl/^An) to make a con- 
nection with the nuclear physics based models, e.g. [10]. The precise value of this 
coupling is not important to us. We know however that it has to vanish when Tq ap- 
proaches r^. This is because the case tq = is a special antipodal configuration, for 
which the world-line of D7— D7 on S*^ degenerates to two antipodal points. In terms 
of parametrization (p = (f){r), the Z2 inside the emergent U{1) would just fiip the sign 
of the fiuctuation 6(f). Therefore in the leading order 6(f) cannot couple to the baryon 
and Qa must vanish. This is not obvious from (5.7), because the expressions we used 
are for the KW case, which is not a good approximation for small tq. Our justification 
of using them to calculate the potential (5.9) in the regime of small ro is that when ro 
is small but not very close to the KW approximation would give a correct answer 
up to an overall factor of order one. To see how the coupling constant vanishes in 
the antipodal case we switch to the KS coordinate t(0), where = e^coshr, and 
remind the reader the form of the induced metric near the tip r = [17]: 



The functions h, Qrr and ga expand as Cq + Cir^ at small r, while g^^ ~ r^. Then the 
coefficient in the interaction term (5.8) becomes proportional to 



where ( is some numerical constant and 6{(f>) is the delta-function. The geometry 
of almost antipodal configuration can be described as follows: r as a function of (f> 
changes very rapidly from infinity to in a very short region near (f> = — 7r/2; then 
it stays constant until (f> almost reaches 7r/2 and after rapidly goes to infinity. This 
means that r is small and constant for almost all values of 0. Near the tip — 0, r 
approaches zero, but r' approaches zero much faster. Therefore expanding (5.11) up 
to the linear order in 6t would not result in any non-zero coupling. 

Similarly we calculate the potential induced by the exchange of the vector mesons. 
Action for the electric part of the U{1) field (3.11) on the D7— D7 branes reads 



3 




(5.10) 




(5.11) 




), (5.12) 
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where the dimensionless Ao(x'^) was expanded in terms of dimensionless functions 
ao(x^) and V9(0). The interaction term comes from the Chern-Simons part of the 
action 

5cs ~ M / d^a;d0 io V {F^JFk^) d^a;aoV^(0) , (5.13) 



with ip{0) ~ 1. The resulting potential is simply 

2 „-m^\x\ 

Kector = /^^^. (5.14) 

AttQs \x\ 

So far we have ignored the pions, the Goldstone bosons of the broken chiral sym- 
metry. These particles are massless and, provided their coupling to baryons is not 
suppressed, the one-pion exchange (OPE) would give a dominant contribution. Thus 
the processes involving OPE are important in the case of interactions of individual 
baryons. The effective pion-nucleon vertex has the form 

2zgA^N{p'h,Yr''N{p). (5.15) 

In the case of massless pions the resulting potential is 

Here Tj is the operator of isospin and the dependence on the operator of spin cTj enters 

Ol2 = 3 5 (Ti ■ (T2 . (5.17) 

A similar potential appears in the Sakai-Sugimoto model [2]. In that case the 
low-energy baryon-pion interaction can also be described by the coupling (5.15) with 
appropriate Qa and The pions vr and massless isoscalar meson t]' come from the 
fluctuations of the gauge field 

A, = U^d,UM<P) , U = exp (^M±ll^ , = (y^^^*)-i . (5.18) 

However the expansion of the D7-brane action around the classical solution for the 
static baryon does not lead to an interaction term. To obtain (5.15) and calculate 
gA one has to quantize the baryons [32]. The resulting potential, besides (5.16), also 
includes an extra term due to the exchange of rj'. The latter coincides with (5.16) 
modulo the isospin interaction term fl ■ T2 and the value of qa- 
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We expect our model to be similar to the Sakai-Sugimoto model. Namely we 
expect the pion (and t]') potential to be of the form (5.16) with appropriate values 
of coefficients and isospin structure. Although we do not compute qa and /^r in this 
work, we expect them to satisfy 

^ — (5 19) 

Clearly the massless pion tail will dominate the large distance potential. At the 
characteristic scale of the a and u potentials, that is x ~ "^gb > the OPE potential 
will be of the same order as the a and u contributions. 

It is believed however that for heavy nuclei or baryonic matter the contribution of 
the OPE is largely averaged to zero due to the tensor structure. Following the nuclear 
physics models of such systems we will ignore its effect on the potential. Since the two 
(or more) particle exchanges correspond to loop diagrams and therefore Qs suppressed, 
we are ready to proceed to a discussion of our findings: in the next section we will 
present a holographic model of weakly bound baryons. 



6 Conclusions. A holographic model of baryons 

In this paper we considered a holographic model of baryons based on non-supersym- 
metric D7 and D7-branes embedded in the KS geometry. We addressed the question 
of baryon interactions and in particular searched for a model of baryons that attract 
at large distances and repel at small. A more specific goal was to shed some light on 
the problem of smallness of the binding energy. These are the identifying properties 
of baryons. 

The holographic approach gave us control over the meson dynamics and made 
other contributions subleading. We have observed that the holographic limit A — 5> oo, 
(7s — > suppressed the glueball contribution to the binding energy, as well as con- 
tribution of the multi-particle exchange. We have also discarded the contribution 
of massless pions by assuming the case of heavy nuclei, or more generally, baryonic 
matter. This leaves us with the lightest massive meson exchanges as the leading in- 
teractions of baryons. The two competing forces are the attraction induced by the 
meson a which corresponds to fiuctuations of the embedding profile and the repulsion 
due to the meson u described by the fiuctuation of the U{1) gauge field on the world 
volume of the D7— D7 branes. 
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In the model we considered, the masses of a and uj are parameterized by rg, the 
position of the lowest point of the D7— D7 embedding. For values of tq far from the 
IR end of the geometry r^, mo- <^ m^^. This ensures the attraction of baryons at 
large distances. The analysis of this paper does not allow to conclude whether the 
attraction would switch for repulsion at small distances, but even if it does so there is 
no natural suppression for the binding energy except for the small factor specific 
to holographic models. 

A more interesting regime is tq ~ r^, when mo- ~ m^. This paper does not give 
the answer, whether mo < always holds, but one can certainly tune the parameter 
ro to satisfy this condition. This sounds promising because m„ ~ m^^ itself implies 
the near-cancelation of the forces. The large distance potential has a typical form for 
effective models of nuclear interactions, e.g. [10]: 



For mo < m^ and < Qw the potential is attractive at large and repulsive at small 
distances. As we discuss in detail in the text, in the regime ro ~ the coupling 
constant (^o- should be small. This is due to an emergent Z2 symmetry, which forces 
g„ to vanish in the case = r^. Therefore one can tune tq to satisfy Qa < Qui- A 
crucial question which needs further investigation is whether one can tune ro to satisfy 
both < Qui and m„ < simultaneously. 

Assuming the answer is positive the resulting model will share the basic property 
of the nucleon physics. The binding energy of two baryons will be suppressed by the 
similarity of the masses of the mesons and potentially by the smallness of the ratio 



Is this a quantitative model of the baryonic interaction? First of all, one has to 
match the mesons a and u to the real world mesons. This problem and in particular, 
the fate of a was extensively discussed in [3]. Next, one can try to match the masses 
and couplings in the potential (6.1) to the real world values. In the present model 
all parameters, including mass ratio and the value of couplings depend on a single 
parameter tq. It seems there is no reason for the matching. Finally, one will have to 
say something about the contribution of pions, multi-particle exchanges and glueballs, 
once the holographic limit is relaxed and they are no longer suppressed. 

To summarize, we formulate the two key ingredients that give us a hope to build 
holographic models of baryon interactions. First, for attractive baryons, one has to 
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satisfy rrifj < m^^. In our case we relied on the conformal symmetry in the UV. 
The meson a is a pseudo-Goldstone boson and hence it is natural that it is light. 
Second, one has to satisfy the condition < g^. In both the Sakai-Sugimoto and 
the present model this can be achieved by dialing tq. At smallest value of tq there 
is an emergent Z2 symmetry which forces the coupling to vanish. The complication 
may arise while combining these two ingredients. The light pseudo-Goldstone boson 
requires the breaking scale tq to be much larger than the internal scale of the gauge 
theory r^. At the same time the Z2 symmetry emerges when tq = r^. The next crucial 
step would be to calculate the masses of a and u for a general configuration to ensure 
nia < rriu]- We leave this task for a future work. 
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A Stability of the D7— D7 embedding 

To show local stability of some classical solution one has to show that the spectrum 
of small perturbations is positively defined. Although simple conceptually, this task 
could be complicated in practice. In the case of interest the main problem would 
come from the interaction of fluctuations of geometrical shape with the world- volume 
gauge fields. Instead we would rather prove stability by showing that the energy 
(action for the static configuration in question) reaches its minimum on the solution 
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of interest. It was shown in [17] that the action of the D7-brane is cahbrated by the 
geometrical volume it spans 



Therefore we are going to show that the volume can only increase should there be 
a small fluctuation around the original solution. This logic allows us to disregard 
the dynamics of the world-volume gauge fields and from now on we simply put them 
to zero. Let us point out that while the world-volume gauge field does not affect 
stability, the effect on stability of the fields in the Minkowski directions is far from 
being trivial [33]. 

In general our goal should be stability of the D7— D7 configuration in the KS case. 
However the corresponding calculation would be quite complicated. Instead we will 
focus on the KW case. When tq ^ the geometry of the KS solution approaches the 
one of the KW solution. We will find that except for the massless Goldstone boson 
a all other modes governing the shape are massive. It is argued in section 2 that the 
Goldstone boson acquires some positive mass when we turn on a non-zero e. This at 
least proves stability in the regime tq ^ r^. 

When ro ~ the KW analysis is not valid. But unless the D7-branes tend to 
spontaneously break the SU{2) symmetry of the KS solution stability follows from the 
uniqueness of the SU{2) invariant classical solution found in [17] and the observation 
in section 2 that the D7— D7 configuration tends to stretch to a smaller radius rather 
then be pulled to infinite r. Although absent in the tq ^ regime, the SU (2)-charged 
instability could theoretically develop when tq ~ r^, but without any evidence for such 
a behavior we consider this possibility to be unlikely. 

In what follows we focus on the stability of the D7— D7 brane in the KW case. 
We start with the singular conifold geometry 




(A.l) 




(A.2) 



where 



= dijj + cos 6d(p + cos 6 d(j) . 



(A.3) 



and consider a static configuration 



e 



(A.4) 



27 



The classical configuration found in [5] corresponds to 



"=2- '°'[7EV=f^- <*-^> 

We want to show that this configuration minimizes (A.l) at least locally. 

A comment is in order. Instead of using the radial coordinate r as an independent 
variable, we are using 0. This is because of the singularity of the r coordinate at 
the tip of the D7— D7 profile. Use of helps to avoid this problem, since the point 
= is regular. The price we pay is mild: the boundary conditions for the SU{2)- 
invariant mode of r are such that it diverges at the end points 0o = ±v^7r/8 as 
r ~ (0 — 0o)~"'^/^, but this will not cause us much trouble. The 5'f/(2)-invariant 
modes were already analyzed in [5] using a different coordinate system. Our goal 
is to focus on the S'?7(2)-charged modes of r and 6, which vanish at the endpoints 
= ±V67r/8. 

For obvious convenience reasons we introduce a G [— 7r/2, 7r/2] through = ^/Qa/A 
and expand the profiles 6 and r as follows 

~9 = ^ + Q{a,ij,9,^), (A.6) 
' R{a,ij,9,^). (A.7) 



cos^/"' a 



The boundary conditions for G and R are such that they vanish at the endpoints 
a = ±7r/2 except for the mode of R, which is independent of ip, 9 and ip. Plugging 
this into (A. 2) and expanding (A.l) up to the quadratic order we get 

5eff = J da dtp d9 d^ ^ ^osingcos^^a ^^^^^,2 ^^^2 ^ ^ 9^2 ^^^2 ^ _ 5^2 

rnsin6'/8^,n 2rn-\/6sinasin6' 



144 V 3 / 144cos3/4a 5^ 

24 sm 9 144 sm 9 cos^ a 



Here we introduced the notations 



(v/f = C = {^,o,i^}, (A.9) 



9'' 



( 1 -cos^ \ 

sin^ 9 
y - cos ^ 1 + 1 sin^ 9 ) 



(A.IO) 
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Let us first say few more tilings about (A. 8) in tfie case wlien tlie perturbations 
R and B liave no dependence on the angular variables ip, 9 and 0. The variation of 
the classical solution (A. 5) with respect to tq gives the zero mode R = (cosa)~^/^ of 
the corresponding equation 

-16^ ( cos^/2^^^ + 9cos'/'' a R- 5cos'/^ aR = Xcos'/^ a R. (A.ll) 
da \ da y 

Clearly, presence of the zero mode implies that the first line of (A. 8) is not positively 
defined. Obviously the problematic term is — 5 cos^/^ ai?^. Later we will show that 
the angular-dependent term cures this problem. In the angular-independent case 
the absence of the negative modes can be derived from the relation between the 
number of zeros and the principal quantum number of the wavefunction of the one- 
dimensional quantum-mechanical problem (A.ll). We would also like to note here 
that the equation (A.ll) is a Heun equation, which admits an extensive symmetry 
group. 

The effective action for the 5'f/(2)-independent modes of 0, 

(A.12) 

is obviously positively defined on the space of functions vanishing at a = ±7r/2. There 
are also zero modes 

e = sin((/) + A(/)) = sin 1^-^ a + A(/)^ , (A.13) 

which break boundary conditions. These zero modes correspond to the rotation of 
S'^ by the broken SU{2) symmetry of the conifold. The constant A0 depends on a 
choice of the axis of the U{1) C SU(2) rotation of S*^. 

Now we introduce the angular dependence and would like to show that (A. 8) is 
positively defined. First we would like to estimate from below the contribution of 
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l [d^dij^ (q cos^/^ aiWRf + {VqA . (A.14) 

14 y sm 6^ V cos2 ' J ^ ' 



First we expand R and B in angular harmonics (recall that ip has periodicity 27r and 
ijj has periodicity 47r): 

= ^ E ^n.,n,(«,^)e^("^'^+"-^) , Rn,,n, = , (A.15) 
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and similarly for G. Here and are either both integer or half-integer. To 
estimate (A. 14) we need to find the ground state of the corresponding eigenvalue 
problem 

' f,inS'^) + <:!i±<Z^hf!^*^^ = X^. (A.16) 



sin ^ d^ V y sin^ 9 

Assuming n^, > the corresponding ground state is 

^ = (cos ^)"^+"^ (sin , (A.17) 

and the corresponding eigenvalue is A = + n, where n = max{n^,n^}. Plugging 
this back we obtain 

167r^rn 
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^ y ^ da cos^/^ a (lQR'i/2,1/2 cos^ a + 9Rl^ ,^^ cos^ a - 5i?^^^„^) + 



^2 _ B'2 - 92 

I 2 n^,n^ ^ntp,n^ 



where 



( 6cos^/^a|i?„^,„^P H ^-^|0n^,n^,M (n^ + 2n - n^i:cos7sina) , (A. 18) 

Y COS Ol I 



cos 7 — 



cos3/4a(6cosV2a|i?„^_^J2 + _lL_|e„^_^j2j ' |K^,„,e*^,„J 

Obviously t, which can always be chosen to be positive, is never larger than 1. The 
term with t in (A. 18) is negative. It causes the largest possible negative contribution 
when C0S7 = 1 and n = n^, i.e. > n^. Since tsina < 1 the effective action (A. 18) 
is positive whenever —5 + 6(n^ + n) > 0. Hence there is only one tricky case to 
analyze: n = n^, = = 1/2 and cos 7 = 1. 

Let us first analyze this case assuming that the perturbation of 6 is zero, that is 
t = 0. Then we have — 5 + 6 ■ (^ + 2|) > and the effective action is positive. What 
happens when t ^ 07 In this case, naively, the action for R is not positively defined. 
If R' is sufficiently small near the endpoints a = ±7r/2, one can expect a negative 
contribution. However there is a positive contribution coming from 0-dependent 
terms that overpower the negative contribution from R. To make the analysis more 
quantitative we estimate (A. 18) from below as follows. First we drop the non-negative 
term R'^ and substitute for (A. 12) its minimal value 50^/3. Then we express in 
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terms of R and t. Eventually we obtain 



n^,n^ 2 



+12 cos^/2 ^ (i ~ ^ ^ ) • 

This is positive for all < t < 1 and — 7r/2 < a < 7r/2, which completes the proof. 

B Size of the baryon 

In the leading approximation we can calculate the mass of the baryon by evaluating 
the tension of the D3-brane wrapping the minimal near the tip of the conifold. 
Indeed, as we argue in section 3, the D3 will prefer to sit at the IR end of the geometry 
for all values of tq. Thus the leading mass of the baryon is always MA, cf. (3.7). 

Let us discuss corrections to the mass in the regime, when the D3 dissolves in 
the D7— D7 to become a flavor instanton. To find corrections and estimate the size 
of the baryon we should find a static instanton configuration on the D7— D7 . The 
procedure is very similar to the calculation done by Sakai and Sugimoto in the case 
of the D8— D8 branes in the type IIA theory [29]. The effective theory for the U{Nf) 
gauge fields on the D7— D7 is given by the DBI part (3.13) and the Chern-Simons 
part (3.9). We split the gauge fields in the f/(l) and SU{Nf) parts (3.11). The non- 
abelian gauge field A will form an instanton configuration charged under the "electric" 
field Aq. At first we ignore the Chern-Simons term, although it is non-trivial due to 
M units of the RR flux through S^, and set the U{1) gauge field to zero A = 0. 

The action for the non-abelian fields A follows from the expression (3.13): 



^ J d'xdz kiz)TT [hFl + 2Fl] , (B.l) 



where k{z) is given by (3.15) and the a radial coordinate z was introduced according 
to (3.14). In terms of this coordinate the geometric profile is z ^ —z symmetric and 
has a parabolic behavior r{z) = rg + 0{z'^) near the tip. 

We are interested in the static (xo-independent) instanton solution with F^i = 
Fqz = 0. Employing the usual trick we find that the mass of such a state is not 
smaller than 



aryon 



> / d'xdz k{z)h'/\z) \e'^'TTF,,F^ . (B.2) 
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Since k(z)h^^'^{z) ~ r is a monotonously increasing function of z we conclude that 
the energy (B.l) is minimal for an infinitely small instanton localized at 2; = 0. The 
mass of such a configuration exactly matches the mass of the D3-brane wrapping the 
minimal (3.7). For an instanton of the small size p, such that p <^ [r^ = e^^^) 
the leading correction to the mass will be of the form 

'''''baryon " '''baryon ^4/3 ' \^-^) 

where a is some positive coefficient of order 1. 

Similarly to the logic of [29] so far we ignored the U{1) gauge field and possible 
contributions from the Chern-Simons term (3.9): 

The instanton carries the U{1) charge and sources the electrostatic field. To calculate 
the corresponding contribution to the baryon mass we will follow the same procedure 
as in the D8— D8 case. We assume that the size of the instanton p is very small 
compared with and therefore we can expand the action around z = effectively 
reducing the geometry to the fiat space. It is also convenient to rescale the spatial 
coordinates Xi Ui = h~^/'^{r^)xi bringing (B.l) to the canonical form such that 
p now is also the spatial size of the instanton in coordinates Ui. We ignore the 
precise values of the numerical coefficients in what follows. Focusing on the static 
configuration with only Aq ^ 0, we have the following effective action 

M jdx'> j d^y dz (Ci A /i(r,) (d,Al + d^l) + C2A0 e^'" Tr . (B.5) 

The energy of the charge of size p scales like l/p^ leading to yet another correction 
toward baryon's mass 



(^"^balyon " "^balyon /^(^^(Ap)^ ' ^^'^^^ 



Here b is some positive coefficient of order 1. Combining this with the correction due 
to curvature (B.3) we obtain for p <^ 

( p2 ^4/3 \ 

mbaryon = MAM +a^ + 6^J +... . (B.7) 



The energy reaches its minimum for 



g2/3 
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Let us analyze this result. First, the condition p holds for large A, hence our 

assumption that p is small is self-consistent. Similarly to the Sakai-Sugimoto case 
the size of the baryon scales like 1/VX and we cannot neglect higher a' corrections 
as their contribution is of order h~^^'^a'p~'^ ~ 1. That is why we are not interested in 
calculating the explicit values for the coefficients a and b, since they will be corrected 
after the higher order terms are taken into account. Yet, we expect the instanton 
size to remain of the order given by equation (3.17). The instanton cannot shrink 
too much due to the large Coulomb energy. Similarly, it cannot significantly expand 
because in this case the energy will increase with the increase of the warped volume 
of S\ 
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